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This paper brings further insight into the recently published A'^-body description of Debye 
shielding and Landau damping [Escande D F, Elskens Y and Doveil F 2014 Plasma Phys. Control. 
Fusion 57 025017). Its fundamental equation for the electrostatic potential is derived in a simpler 
and more rigorous way. Various physical consequences of the new approach are discussed, and this 
approach is compared with the seminal one by Pines and Bohm [Pines D and Bohm D 1952 Phys. 
Rev. 85 338-353]. 
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I. INTRODUCTION 


In order to discuss microscopic plasma physics, most textbooks start with invoking the 7V-body description of 
plasmas, but they deem it too difficult to tackle, and they rapidly switch to kinetic equations. The kinetic approach has 
proved extremely powerful to compute a wealth of plasma effects, but it often hides some aspects of their mechanical 
nature ; for instance, after the Vlasovian derivation of Landau damping, most textbooks use simple mechanical models 
to provide physical interpretations of this damping. 

However, the A-body description of plasmas has proved for more than six decades to be very powerful at describing 
analytically effects like Debye shielding, Landau damping, spontaneous emission of Langmuir waves, and at revealing 
the corresponding mechanical behavior of both particles and waves. This started in 1952 with a seminal paper by 
Pines and Bohm From the nineties on, there was a revival of the A-body approach to describe analytically 

wave-particle interaction in plasmas [3, [l^ . 

Last year, reference [l^ provided a direct and simultaneous derivation of Debye shielding and Landau damping by 
an A-body description of the plasma. This derivation works directly with Newton’s second law and with a version of 
the electrostatic potential linearized from ballistic orbits of the particles ; it does not use probabilistic arguments or 
partial differential equations. However, in retrospect, the formulation of this paper now appears to be unnecessarily 
intricate, because it was still too close to the heuristic derivation of the theory ; it also required an erratum [l^ . 

The heuristic derivation substituted the discrete summations over the particles, which are intrinsic to the A- body 
approach, with integrals over a continuous velocity distribution, and was presented at the beginning of reference [iSjj . 
In retrospect, this substitution turns out to be an efficient shortcut to the simultaneous derivation of Debye shielding 
and Landau damping, but this non rigorous step can be avoided by considering the initial A-body distribution as 
close to a set of monokinetic arrays of particles. This other approach was presented in a second part of reference 
. In strong resonance with Dawson’s seminal multi-beam fluid theory Q, this second approach shows that Landau 
damping is due to phase mixing a la van Kampen [l5l| . while the first approach leads to the Vlasovian expression 
introduced by Landau that hides this fact. 

Sections El and mil (and the appendix) provide a shorter and more rigorous version of the approach of reference 
using a smooth distribution function. The calculation reduces to a mere exercise using Fourier series and Laplace 
transform. Section |TV] discusses the various physical consequences of the present A-body approach. Section |V] is 
devoted to the comparison of the present A-body approach with the 1952 one by Pines and Bohm This latter 
work is often overlooked by plasma physics textbooks and was overlooked also by the present authors in reference 
M- The main technical difference of the present theory with reference is the use of the Laplace transform. It 
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turns out to have the same heuristic power as in the Vlasovian approach, and enables a compact and straightforward 
simultaneous derivation of Debye shielding and Landau damping without any a priori physical intuition. 


II. SIMPLER DERIVATION OF THE FUNDAMENTAL EQUATION FOR THE ELECTROSTATIC 

POTENTIAL 

As in reference [l^, we deal with the One Component Plasma (OCP) model [U, 0, which considers the plasma 
as infinite with spatial periodicity L in three orthogonal directions with coordinates {x,y^z), and made up of N 
electrons in each elementary cube with volume L^. Ions are present only as a uniform neutralizing background, 
enabling periodic boundary conditions. We now provide a shorter and more rigorous version of the derivation of 
reference using a smooth distribution function. The actual calculation involves sums over the N particles and 
over all Fourier components of the potential. However, the principle of this calculation can be explained more simply 
by performing it formally for the unphysical case of a single electron acted upon by a single Fourier harmonic of its 
own held. This is done in the next subsection, while the complete derivation is in appendix. 


A. Single particle acted upon by a single Fourier harmonic of its own field 


While the Fourier decomposition of the periodic delta distribution corresponding to the charge density of an electron 
entails an inhnity of harmonics, in this subsection we keep just a single one with wavevector k which provides through 
Poisson equation the single harmonic of the electrostatic potential 

‘fit) = exp[-ik ■ r(t)], (1) 

where k = ||k||, cq is the vacuum permittivity, and r(t) is the position at time t of the particle acting as a source with 
charge —e. Let Tq and v be respectively the initial position and velocity of the particle, and let Ar(t) — r(t) — rg — vt. 
We now compute a perturbative solution to the dynamics. To this end, setting r(t) — rg + vt + Ar(t) in equation 
o, we replace (p with its expansion to first order in Ar(t) 

<P\init) = exp[-ik • (ro + vt)] [1 - ik • Ar(t)]. (2) 

We now use the time Laplace transform which maps a function g{t) to 'g{oj) = g{t) exp(ia;t)dt (with w complex). 
Since the arguments of functions are spelled explicitly, from now on we omit diacritics for the Laplace (or Fourier) 
transformed quantities. The time Laplace transform of equation Q is 


‘/5iin(w) =-— exp[-ik ■ ro] [- 

eo/c i 


-k- 


— ik • Ar{uj — k ■ v)]. 


where the Doppler shift —k ■ v comes from the linear dependence on t of the exponent of equation ©■ 
To compute Ar(uj — k • v), we use Newton’s equation, which writes 


( 3 ) 


r = —V(pii„(r) 


TOe 


( 4 ) 


for the particle with mass me- The single harmonic of the electrostatic potential defined by equation © implies 


‘/?iin(i’) = jg <^iin exp(ik • r), 


( 5 ) 


where r now is the position of the particle as a massive object subjected to a force. 

In order to derive our perturbative solution, we consider A:||Ar|| to be small in equation so that the linearized 
particle dynamics defined by equation (|T]) reduces to 


ie 

L^rUf. 


Air = 


k (^iin(t) exp[ik • (ro + vf)]. 


( 6 ) 
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Since Ar(0) = 0 and Ar(0) = 0, the time Laplace transform of equation ([6]) is 


- w^Ar(w) = ^ 3 ^ k exp(ik • Tq) (/3iin(w + k • v), (7) 

where the Doppler shift k • v comes from the linear dependence on t in the exponent of equation Computing 

Ar(a; — k • v) in equation ([3]) from the expression of Ar{uj) given by equation ([7]) yields 


k^(piin{uj) - 


nieCoL^ 


kk 


V?iin(a' + k ■ V - k ■ v) 
(w — k • v)2 


exp[i(k - k) • ro] = k'^, 


( 8 ) 


where 


(bal) / \ 


ie exp[—ik • rg] 
eok^ w — k ■ V 


(9) 


is the time Laplace transform of the ballistic potential obtained by setting Ar{t) = 0 in equation ([5]). We keep on 
purpose the terms (k — k) • rg and k • v — k • v due to the two Doppler shifts, in order to exhibit the structure of the 
equation when the full set of wavevectors is taken into account. 


B. Fundamental equation for the electrostatic potential 

As derived in the appendix, when the full set of wavevectors and particles is taken into account, equation (|S]) 
becomes 


w)- 


N 


km-kn 


n, Ain ^smooth ^ 1 


N 


(n,a; + k„ 


■ Vj - K 


■V,) 


exp[i(k„-km)-rjo] = (10) 


-k„ 


1? 


where Wp = [e^n/(meeg)]^/^ is the plasma frequency (with n = N/L^), where y!iin(ni, w) is the time Laplace transform 
of the Fourier component of the electrostatic potential with wavevector km = ^ m, and m = {mx,Tny,mz) is a vector 
with three integer components (we comment below on the restriction fc„6smooth < 1 making kn run over a bounded 
domain). Note that = g(—m, —w*) for real m and complex uj if g{r,t) is real-valued. We also defined 




(bal) / X 


with 


(bal) / X 

7’iin./(m,w) 


ie exp[-ikm ■ fog] 
egfc^ w - km ■ Vj 


( 11 ) 


( 12 ) 


the ballistic potential related to particle j, i.e. the unshielded Coulomb potential of an electron with position rj{t) = 

rjo -f Vjt. 

In this complete derivation, for any finite value of the ||Ar||’s due to the various particles, the linearization in 
the Ar’s i^ustified for finite values of the wavenumber moduli fcm = ||km|| only. Therefore, as in section 4.5 of 
reference [l5|, we restrict the Fourier expansion of the Coulomb potential to fcm’s such that fcm ^smooth £ 1, where 
^smooth ^ Ad = [eofcBT/(ne^)]^/^, the Debye length, with fce the Boltzmann constant, and T the temperature ; let 
us denote by Adgmooth the set of such wavevectors. Then, we may consider finite ||Ar||’s, with ||Ar|| <C ^smooth- The 
above truncation of the Fourier expansion smoothes the Coulomb potential. We notice that a similar smoothing of 
the 1/r singularity of this potential (or some penalization of very close approaches [1^) is performed in the mean-field 
derivation of the Vlasov equation i i [n m M m. In any case, the Vlasov equation cannot describe scales 
smaller than the inter-particle distance. Since the self-field due to the smoothed Coulomb potential vanishes, it is not 
necessary to exclude self-interaction in the complete derivation [l^ . 

Equation (IIOI) is the fundamental equation of this paper. This equation is of the type £ip\i-a = source term, where £ 
is a linear operator, acting on the infinite dimensional array whose components are all the Doppler shifted (piin(ni, w)’s. 
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III. SMOOTHED LINEAR OPERATOR OF THE POTENTIAL EQUATION 


We now approximate the granular (or empirical, Kolmogorov) distribution F{r,v) = — r^o) (5(v — v^) 

with a position and velocity distribution function /(r, v) that is continuous in r, such that distribution f — F yields a 
negligible contribution when applied to space dependent functions which vary slowly on the scale fosmooth- This means 
that we use a version of F coarse-grained on this scale. On replacing the discrete sums over particles with integrals 
over the smooth distribution function /(r,v) in its left hand side, equation (flOl) becomes 


w) - ^ ^ • k, 


(/;cg(n, cu -I- (kn - km) ■ v) 

(w - km • v )2 


/(k„ -km,v! 


d^v = 


(13) 


where tpcg is the smoothed version of </5iin, and / is the spatial Fourier transform of /. We consider almost uniform 
initial spatial distributions of particles. So, / is a spatially uniform distribution function /o(v) (fj /o(v)d^rd^v = 
f /o(v) d^v = L^n = N) plus a small perturbation /i(r, v) of the order of if eg- In linearizing equation (IT^ for fcg, 
/(k, v) is substituted with /o(k,v) which vanishes for k 7 ^ 0. Therefore, operator £ becomes diagonal with respect 
to both m and uj (a complex quantity), and equation (fT51) becomes 


with 


e(m,w)(/Jcg(m,w) = (/3[i'(f'^(m,a;), 


, i2r3 

e(m,a;) = 1 — ^ 


/o(v) 




(14) 


(15) 


The validity of the above coarse-graining was discussed in section 4.5 of reference [T^. Here we only recall the 
main steps of the discussion. First, the presence of many particles in the Debye sphere justifies neglecting the non¬ 
diagonal terms in equation OT. Second, Dawson’s work Q justifies the Vlasovian expression m for e(m, w) for a 
one-dimensional plasma. The generalization of Dawson’s calculation to the three-dimensional case, though tedious, 
does not a priori involve any new conceptual difficulty. These results back up the above coarse-graining. 

However, we now show these results come with conditions on the parameters of our problem. Invoking Dawson’s 
calculation means that, when interested in a wave with phase velocity (with ||vw|| = v^), there are many particles 
with parallel velocities in the range [uw — Iqul/fc, ^’w + Ibul/fc], where 7 l is the Landau growth or damping rate of the 
wave, and k the modulus of its wavenumber. Indeed, this range is the one where the phase mixing a la van Kampen 
is occurring, and where the synchronization of particles with the wave brings the change of particle momentum 
inducing Landau damping or growth of this wave I,EE in,[13 ■ If the initial distribution of particles is meant to 
approximate a smooth distribution /o(v), this means that the length L of the system is taken large enough to fulfill 
the condition “many particles with parallel velocities in the range [ww — ItlI/^j'^w + ItlI/^]”- However, whatever L, 
for V large enough, one should not use e(m,a;), but its discrete counterpart computed with the A-body distribution 
function. 


IV. DEBYE SHIELDING, LANGMUIR WAVES AND LANDAU DAMPING 

To keep this paper self-contained, we now discuss how Debye shielding and Landau damping were derived from 
equations (fTS - dTSl) in reference [l3. These equations are obtained by introducing the continuous velocity distribution 
after particle dynamics is taken into account, and not before, as occurs when kinetic equations are used. However, 
they look very similar to those including initial conditions in Landau contour calculations of Langmuir wave growth or 
damping, usually obtained by linearizing Vlasov equation and using Fourier-Laplace transform, as described in many 
textbooks. Indeed, the smoothed self-consistent potential fcg is determined by the response function e(m,a;), viz. the 
classical plasma dielectric function, whose Vlasovian expression involving dfo/dv obtains by a mere integration by 
parts if /o is differentiable. For a cold plasma, e(m, w) = 1 — Wp/w^, as expected (however this case is ruled out since 
the present theory is correct only for a large number of particles in the Debye sphere, as recorded in the previous 
section). 

In equation (fCT) . (m, w) is a discrete sum over the N particles. As computed in [1, EiUll) in the inverse 
Laplace transform of (^cg(ni, w), the pole w = k^ • (0) of (m, w) brings a contribution to (^cg(ni, t) which is the 

shielded Coulomb potential of particle j. This potential is often computed in textbooks by adding a test particle to a 
Vlasovian plasma (see for instance chapter 9 of [13). Since Debye shielding is effective over distances larger than the 
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Debye length, this discussion is relevant under the assumption that our box size is large in comparison, viz. L ^ Ad- 
As in the Vlasovian approach, the part of (/5cg(m, t) due to dominant poles of l/e(m, w) corresponds to one Langmuir 
wave propagating in the direction of and one propagating in the opposite direction, with the corresponding Landau 
growth or damping. Therefore, the full potential in the plasma turns out to be the sum of two contributions only : 
one due to Langmuir waves and one due to all the shielded Coulomb potentials of individual particles. 

Equation (fT4)l becomes the Vlasovian one on replacing in equation (fT^ the discrete sums over particles with integrals 
over the smooth distribution function /(r,v). This yields 




(bal) 

eg 


(m,a;) 


r /i(m,v) 
eok^ J w - km ■ V 


(16) 


the Vlasovian formula for the ballistic potential. The AZ-body description reveals that the contribution of (m, ai) 
to (pcg{m.,t) in file Vlasovian approach is the continuous limit of the sum of all the Debye-shielded potential of the 
particles. In this approach, this contribution is only a transient [^ . Though its physical meaning is unclear, it does 
not correspond to a robust supplementary contribution to the full electrostatic potential in the plasma which would add 
to the potential of Langmuir waves and to the shielded Coulomb potentials of particles. 

To afford the description of nonlinear effects in wave-particle dynamics, section 6.1 of reference [l^ generalizes 
the linear analysis analogous to that in section |n] of the present paper by applying it to the particles which are not 
resonant with Langmuir waves, while keeping the exact contribution to the electrostatic potential of the remaining 
N' tail particles that may experience trapping or chaotic dynamics. This eventually leads to equations describing 
the self-consistent dynamics of M Langmuir waves and of N' tail particles (these equations could be also obtained 
with the formalism of the present paper). For a one-dimensional plasma, the linear theory of Langmuir waves can be 
formulated on this basis, with the bonus that the information on particle dy namics is provided in parallel with the 
waves’ This approach also provides spontaneous emission [12, [H, Il3 ; and some nonlinear effects like 

damping with trapping [l8j |. and the decoupling of the waves from the strongly chaotic motion of the particles when 
there is a plateau in the particle velocity distribution function (see sec. 2.2 of Q)- 


V. COMPARISON WITH PINES AND ROHM 1952 

In 1952, Pines and Bohm showed that an A-body approach can bring a lot of insight into basic plasma effects 
[131. We now summarize the part of their results overlapping with our theory. Pines and Bohm already considered a 
One Component Plasma in a periodic box and wrote the equations of motion of its N electrons acted upon by their 
electrostatic potential. Working with the physical time t, they got an equation (number (8) in their paper) for the 
time evolution of a Fourier component pk of the charge density. This equation involves a contribution of all Fourier 
components whose structure is similar to equation but without linearizing particle dynamics. In order to neglect 
the contribution of all other Fourier components, they introduced the now celebrated “random phase approximation” 
by showing that the particle positions may be considered as nearly random for all times if there are many particles 
in a Debye sphere of a plasma in thermal equilibrium. 

As a first step, they considered the collective part of pk, and showed that, for small wavenumbers ||k||, the charge 
density oscillates at the plasma frequency, and that plasma oscillations exist for ||k|| ^ Inductively, they derived 
a harmonic equation for non-vanishing ||k||’s with an oscillation frequency w assumed implicitly real, and satisfying a 
discrete version of the Vlasovian dispersion relation (their equation (23)). They recovered the Bohm-Gross dispersion 
relation for small ||k||’s. 

As a second step, in their section III they split pk into a collective and a granular part. The evaluation of the latter 
for thermal particles yielded correctly the Yukawa version of Debye shielding. 

In contrast, as those presented in most textbooks, our approach to Debye shielding and Landau damping is a linear 
theory. Laplace transform can then be applied readily. As already demonstrated by Landau’s theory, the Laplace 
transform is a powerful heuristic tool. In the A-body approach, it shows the full potential in the plasma to split 
naturally in two contributions related to two families of poles : one due to Langmuir waves subjected to Landau 
growth and damping, and one due to all the shielded Coulomb potentials of individual particles, whatever be their 
velocity. By considering the initial A-body distribution as close to a set of monokinetic arrays of particles, it was 
shown in Appendix B of reference [l3 | that linearization provides naturally the cancellation of non-diagonal terms in 
equation m ; it is not necessary to use the random phase approximation, or other statistical arguments. 

Particle deflection and Gauss’ theorem were used by Pines and Bohm to interpret their Yukawa version of Debye 
shielding. Section 5 of reference [l5l| brings more mechanical insight to this interpretation. There, instead of applying 
the Laplace transform to linearized particle dynamics, one applies Picard iteration technique (a standard method to 
prove the existence and uniqueness of solutions to hrst-order equations with given initial conditions) to the equation 
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of motion of a particle, say P, due to the Coulomb forces of all other ones without linearizing particle dynamics. 
This calculation shows that particle P acts on another particle P', directly by its Coulombian force, and indirectly 
by deflecting the orbits of all other particles whose force on particle P' is modified by this deflection. The direct 
Coulombian force is reduced by the indirect action, since the Coulombian deflections due to particle P diminish the 
number of other particles in a sphere surrounding it, which diminishes the apparent charge of particle P according to 
Gauss’ theorem. This shielding effect becomes dominant at a distance on the order of Ad- 

When starting from random particle positions, the typical time-scale for shielding to set in is the time for a 
thermal particle to cross a Debye sphere, i.e. 0 ;“^. This is also the time for the equilibrium pair-correlation function 
to settle at its time-asymptotic values. Therefore, shielding results from the accumulation of almost independent 
repulsive deflections with the same qualitative impact on the effective electric field of particle P. If point-like ions 
were present, the attractive deflection of charges with opposite signs would have the same effect. The Coulombian 
deflections induce a self-organization of the plasma : it becomes a dielectric. So, these deflections do not produce 
only the disorder leading to collisional transport, but also the order corresponding to shielding and to the behaviour 
as a dielectric. Paradoxically, “collisionless” Landau damping turns out to occur because of what is usually called 
“collisions”. Since the global deflection of particles includes the contributions of many other ones, the density of the 
electrons does not change, at variance with what occurs for the shielding at work next to a Langmuir probe. Gauss’ 
theorem was used in section 2.2.1 of to provide an intuitive explanation of the shielding of such a probe. 


VI. CONCLUSION 

This paper provided further insight into the unified A^-body description of Debye shielding and Landau damping. 
First, by deriving the fundamental equation for the electrostatic potential of reference in a simpler and more 
rigorous way. Second, by stressing that the term from initial conditions in Landau’s calculation of Langmuir waves 
does not bring new physics, but is a mere damped continuous version of the Debye shielded potentials of all particles 
whose physical meaning is unclear in the Vlasov picture. Third, by discussing more precisely the condition for 
approximating the Wbody distribution by a smooth one in the dispersion relation. Fourth, by comparing the present 
approach with the 1952 one by Pines and Bohm which already brought a lot of insight into Debye shielding and 
Langmuir waves. In particular, we have shown how the Coulombian deflections of particles bring order to the plasma 
by producing Debye shielding and by making the plasma behave as a dielectric. Furthermore, the Laplace transform 
confirms to be a powerful heuristic tool providing a compact and straightforward simultaneous derivation of Debye 
shielding and Landau damping. Linearization was shown to provide naturally the cancellation of non-diagonal terms 
in the fundamental equation for the electrostatic potential. This avoids using the random phase approximation or 
statistical arguments. 

The iV-body approach provides a derivation of Debye shielding for a single mechanical realization of the plasma. 
This was used in an V-body calculation of collisional transport of the One Component Plasma (OOP) incorporating 
all impact parameters with no ad hoc cut-off [l3| . This calculation explains why a two-body calculation yields a correct 
estimate of collisional transport, while most of this transport is due to the simultaneous action of many particles with 
impact parameters between the inter-particle distance and the Debye length. 

When generalizing this theory to the collisional transport of a plasma with granular electrons and ions, the sin¬ 
gularity of head-on collisions of particles with opposite charges can be avoided by excluding the case of a vanishing 
impact parameter which has probability 0. Then, the deflection of particles with same sign is just the opposite of 
the deflection of particles with opposite signs, so that the classical Rutherford cross-section is the same for charges of 
like sign and opposite signs. However, one might wonder whether a quantum mechanical treatment would not modify 
the scattering cross-section. Yet, it does not. This result, traditionally derived with the Born approximation, was 
obtained more recently with a much simpler calculation using Fermi golden rule that makes intuitive why the classical 
and quantum results agree [ 2 ^. 

The theory of collisional transport of reference [ij and the approach of the present paper show that 7V-body 
dynamics, which has always been the ultimate reference in plasma textbooks, is now also a practical tool. Furthermore, 
Laplace’s dream was not a mere utopia, since classical mechanics can genuinely describe non trivial aspects of the 
macroscopic dynamics of a many-body system. 

One of us (DFE) thanks the Program Gommittee of the 42st EPS Conference on Plasma Physics and G.L. Delzanno 
to have led him to look for a more pedagogical derivation of the results of reference [isj] in order to prepare an invited 
talk and a course ; this was the origin of section [IT] of this paper. He also thanks J. Daligault for recording him 
the pioneering 1952 contribution of Pines and Bohm to the Wbody approach to basic plasma physics, and 
J. Daligault and S. Khrapak for a critical reading of our manuscript. We thank the two anonymous referees for 
leading us to formulate our theory in a more rigorous way, and to write a clearer and more self-contained paper. 
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Appendix A: Full derivation 


We now consider (p{r), the potential created by the N particles at any point where there is no particle. Its discrete 
Fourier transform is readily obtained from the Poisson equation, and is given by (^(0) = 0, and for m 7 ^ 0 by 


^(m) 


eokl 


N 

■^exp[-ik 

r 

i=i 




(Al) 


where rj{t) is the position at time t of particle j acting as a source, (^(m) = J (p{r) exp(—ikm -r) d^r, with km = ^ Hi) 
and km — ||km||- 

Let Tjo and Vj be respectively the initial position and velocity of particle j, and let Arj{t) — rj{t) — rjo — Vjt. 
We now compute a perturbative solution to the full fV-body dynamics and the resulting Fourier components of the 
potential for m running in the finite domain Alsmooth defined in subsection III B[ such that fcm&smooth < 1- 

To this end, setting rj = Vjo + Vjt + Arj{t) in equation (IA1|) . we replace (p with its expansion to first order in the 

Arjitys 


N 

(^iin(m,t) = —p- exp[-ikm ■ (r^o + Vjt)] [1 - ik^ ■ Arj{t)], 


j=i 


whose time Laplace transform is 


N 


99iin(m,w) = - W —^ exp[-ikm • r^o] [-r-ikm • Arj{u} - km • v^)], 

7=1 ^O^m ^ ~ ■ Vj 


(A2) 


(A3) 


where the Doppler shift —km ■ Vj comes from the linear dependence on t of the exponent of equation (IA2I) . 

To compute Arj{uj — km ■ Vj), we use Newton’s equation for the particles 

icj = ^Vpiinirj), (A4) 

TiLq 

where rrie is the electron mass. When m runs over the finite domain such that fcm&smooth £ 1) the self-field due to (/j 
vanishes, and it is not necessary to exclude self-interactions. Therefore, one may use the harmonics of the electrostatic 
potential due to all particles defined by equation ini, which yields 


(^ii„(m)exp(ikm • r^). 


(A5) 


where rj now is the position of particle j as a massive object subjected to a force. Using equation (IA2I1 . the linearized 
particles dynamics defined by equation (Ell is then given by 


Arj = j^— kn (^iin(n,t)exp[ik„ • (rjo + Vjt)], (A 6 ) 

n, ^‘n^smooth ^ 1 

Since Arj = 0 and Arj = 0 at t = 0, the time Laplace transform of equation (IA 6 p is 

- uj‘^Arj{uj) = y^ k„exp(ik„ • Tjo) </ 5 iin(n,a; + k„ • Vj), (A7) 

n, fcnf*smooth ^ 1 

where the Doppler shift kn • Vj comes from the linear dependence on t in the exponent of equation (IA 6 II . Computing 
Ar{u! — k • v) in equation (IA3() from the expression of Ar{uj) given by equation JAZI yields the fundamental equation 
(fTUl) for the potential. Note that the assumption fcm&smooth < with ^smooth Ad, defining the set Adsmooth, excludes 
scales which are irrelevant to Debye shielding and Landau damping, since these phenomena involve scales larger than, 










or of the order of Ad • 
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